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Abstract 

We reanalyze the question of black hole creation in high energy 
scattering via Shockwave collisions. We find that string corrections 
tend to increase the scattering cross-section. We analyze correc- 
tions in a more physical setting, of Randall-Sundrum type and of 
higher dimensionality. We also analyze the scattering inside AdS 
backgrounds. 
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1 Introduction 



The problem of black hole creation in high energy scattering is one of significant importance, 
for two possible reasons. One is that one can have a low gravitational scale, as in the 
large extra dimensions 12] or Randall-Sundrum [2 0] scenarios. Thus the possibility of 
black hole creation at accelerators has been explored at length in a number of papers (e.g. 

ia in ei ca im ca C2i ) • 

Another reason, is that via gauge-gravity dualities, the high energy scattering in a gravity 
theory can be related to high energy scattering in QCD, or gauge theories in general [T4* l ll5 | 
ITS] . Simply put, high energy scattering in QCD can be described in terms of a conformal 
field theory with a cutoff, and that is dual to a two brane Randall-Sundrum scenario. But 
then black hole creation that happens in high energy gravity scattering has to have some 
implications for the QCD side. In fact, in ^Hj it was argued that black hole creation, when 
the black hole size is comparable to the size of the gravity dual= AdS slice, is responsible for 
the much sought-for Froissart behaviour (saturation of the unitarity bound). We will revisit 
these questions in a future paper [T7j, but we will still set up some of the calculations needed 
for that case in here. 

In particular, we will analyze the case of high energy scattering with black hole formation 
inside AdS space. 

We will focus instead on the actual black hole creation at high energy s ~ Mp t with the 
idea of applying it to theories with a low fundamental scale. 

Giddings and Thomas [3] and a number of other people jSJ El E3 GDI H2j (see also earlier 
work in jjj) have proposed that the cross-section for black hole creation in flat space at high 
energy is just proportional to the geometric horizon area of a black hole of mass equal to the 
total center of mass energy, i.e. 

a~nr 2 H ; r H = 2G^~s (D = 4) (1.1) 

There has been a considerable amount of debate over whether this assumption is correct 
(see, e.g. EH CHI CH C3 EH ) • 

In an attempt to prove it, Eardley and Giddings have treated the high energy collision 
according to a recipe proposed some time ago by 't Hooft [20] • The process is well described 
by the collision of two gravitational Shockwaves of Aichelburg-Sexl type. Even though one 
cannot calculate precisely the metric in the future of the collision except perturbatively |21j . 
one can use a trick due to Penrose that just uses the properties of Einstein gravity to calculate 
a lower bound on the area of the horizon that will form in the collision. 

In D=4 were able to extend Penrose's method to collision at nonzero impact pa- 
rameter b of the two Aichelburg-Sexl waves, and prove that the cross section for black hole 
scattering is indeed of the order of magnitude of the geometric cross-section of the classical 
black hole. 

In this paper we will try to refine this calculation, and answer some of the criticisms 
addressed to the calculation and the geometric cross section result. One such criticism was 
that string corrections will significantly lower this result (see [22] for example) We will try 
to analyze string corrections explicity via two methods. 



2 



There are two modifications of the Aichelburg-Sexl metric that were shown to reproduce 
string scattering results (effective metrics). The one in [22] analyzes specifically the scattering 
at impact parameter b, and gives an effective metric for large b (> R s , the gravitational 
radius for black hole formation). It is therefore unsuited for our purposes, yet with some 
approximations one can find that the head-on collison of two such waves (each having a 
parameter b) will have an increased horizon area of the formed black hole, with respect to 
the Aichelburg-Sexl case. The second modification [21] corresponds to string-corrected 't 
Hooft scattering in an Aichelburg-Sexl metric. We will show that scattering of two modified 
Shockwaves will again increase the horizon area of the formed black hole. 

Another possible caveat to the calculation in ^T] is that it was done in flat D=4. We will 
analyze the case of the more realistic Randall-Sundrum scenario and find that we just get 
small corrections to the flat D=4 case. We will also offer a method of estimating the cross 
section in the arbitrary D case. 

The paper is organized as follows. In section 2 we will review the Aichelburg-Sexl wave 
and 't Hooft's scattering calculation, and generalize it to higher dimension. In section 3 
we will review the analysis of |TT] and set it up for generalization to any Shockwaves and 
any dimension. We will also analyze the collision of sourceless waves, which should describe 
graviton-graviton scattering, and present a puzzle. In section 4, we will analyze string 
corrections via the effective metrics in [23] and [21]- In section 5 we analyze the case of 
Randall-Sundrum background and calculate corrections. In section 6 we will write down a 
solution for an Aichelburg-Sexl wave inside AdS and do a 't Hooft scattering analysis. 

2 The Aichelburg-Sexl wave and 't Hooft scattering at 
high energy 

't Hooft [20] has proposed that an (almost) massless particle at high energies s ~ Mp t ^> t 
behaves like a plane gravitational wave- a shockwave, and its only interactions are given by 
massless particles, with the gravitational interactions described by deflection in the gravita- 
tional shockwave corresponding to the massless particle. That shockwave solution is due to 
Aichelburg and Sexl |25j . 

In this section we will review this procedure of gravitational interaction and generalize 
it to higher dimensions. 

The Aichelburg-Sexl solution is of the pp wave type. A pp wave (plane fronted gravita- 
tional waves) has the general form in d dimensions 



and the rest are zero. Horowitz and Steif j2H] showed that there are no quantum (a') 
corrections to the (purely gravitational and NS-NS background) pp wave solutions, since all 



ds 2 = -dx + dx~ + (dx + ) 2 H(x + , x l ) + ^(dxi) 2 




and has Ricci tensor 
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the gravitational invariants made from Ricci and Riemann tensors vanish on this solution. 
The inverse metric is given by 

g^dypv = -4<9+<9_ - 4Hd 2 _ + d 2 (2.3) 

and so for instance 

= R^BT = R^R^g™ (2.4) 

does not contain (R ++ ) 2 , and is thus zero. 

In particular, a class of purely gravitational (sourceless) solutions (of R ++ = 0) are given 

by 

H = Y^ AjX l x j , trA = (2.5) 

ij 

and preserve 1/2 susy T_e = 0. 

The Aichelburg-Sexl solution is a solution for a point particle (delta function source), 
moving at the speed of light. It is obtained by boosting the black hole solution to the 
speed of light, and taking its mass M to zero, while keeping Me@ = p = const. (/9=boost 
parameter). But a simpler way to get it is to boost the energy momentum tensor and then 
solve the Einstein equations for the resulting pp wave (thus we have to assume the pp wave 
ansatz, which however turns out to be consistent with the energy- momentum tensor). 

A black hole at rest has 

T 00 = m Q 5 d - 2 (x l )5(y) (2.6) 
and the rest zero. Boosted, one gets 

T <n = -jp=5*-*(z i )5( y -vt) (2.7) 
V 1 — u 

and corresponding T 10 and T n . At the limit, one has 

T ++ = p5 d ~ 2 (x i )5(x + ) (2.8) 

This means that H(x + ,x l ) = 5(x + )Q(x' 1 ), where (since Einstein's equation is -R++ = 
8nGT ++ ) 

d 2 ^(x l ) = -l&nGpS^tf) (2.9) 

($ is harmonic with source). 

For 4d gravity, $ = SpGlnp, and 

ds 2 = -dudv - ApGlnp 2 5(u)du 2 + dx 2 + dy 2 (2.10) 

in the notation of |2Z| {p 2 = x 2 + y 2 ), but the result is easily generalizable to any dimension 
d higher than four: 

$=- -r-f-r, d>A 2.11 

Particles following geodesies in the A-S metric are subject to two effects [27] : 
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It is found that geodesies going along u at fixed v are straight except at u=0 where there 
is a discontinuity 

P 2 

Av = $ = -AGp ln£- (2.12) 

where the Planck constant Ipi in the In is conventional (only relative shifts, Av\ — At>2, 
have physical meaning). That means that one basically has two portions of flat space glued 
together along u=0 with a Av shift. The shift can be easily understood by the fact that after 
a singular coordinate transformation, defined later on in (|3.4j) . the metric becomes continous. 
So geodesies are continous in (u, v) coordinates, which means they are discontinous in (u, v) 
coordinates, with the above Av. 

The second effect is a "refraction" (or gravitational deflection, rather), where the angles 
a and j3 made by the incoming and outgoing waves with the plane u = at an impact 
parameter p = b from the origin in transverse space satisfy 

AGv 

cota + cot (3 = (2.13) 

(here p is the momentum of the photon creating the A-S wave), and at small deflection 
angles (near normal to the plane of the wave) we have 

AO ~ — ^ (2.14) 
b 

We can understand this also by using the singular coordinate transformation in ([3.4)1 . as 

<"»> 

and A(|£) = (no refraction in (p,u,v) coordinates), so 

One can then describe the scattering of two massless particles of very high energy [20] 
(mi,2 Mp, Gs ~ 1, yet Gs < 1) by saying that particle two creates a massless shockwave 

(2) 

of momentum and particle one follows a massless geodesic in that metric. In covariant 
notation (v = x° — x 1 = x~, x 2 = p 2 = x 2 + y 2 ), 

A Xfl = -2Gp {2) log(x 2 /C) (2.17) 

Then particle one comes in with a free wavefunction 

,0(1) _ e ipWx+ip"'v+ip^u ^2 18) 

and becomes (at u=0, just after the shockwave) 

= e ip m S=+ip ( -\v-iGp^log{x 2 /C)) {2.19) 
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Then by definition the scattering amplitude is the Fourier transform of this wavefunction 



(27T) 



= -i 5 (k^ -pW) [ J^Le^e^ (2.20) 
J (2tiT 

where we have expressed A(s, t) via an impact parameter transform to an eikonal form, with 
5(b, s) = p+?Av = —Gslogb 2 and after doing the d 2 b = bdbdO integration one gets 't Hooft's 
result 

But 

AGp [ l ] p {2) = Gs and (p - k) 2 = -t; A{k + , k) = 5{k { + } - p { + ] )U(s, t) (2.22) 
and then we get the differential cross-section 

UM = i ( ± )M0 .r(i-«J.) 



4n y -t' T(iGs) 

which is like Rutherford scattering, as if a single graviton is exchanged, (with the effective 
gravitational coupling Gs replacing a = e 2 /4ir of QED) 

The argument is that graviton exchange dominates the amplitude in this limit, for massive 
particles it takes an infinite time to interact. Indeed, at large impact parameter there is the 
natural exponential decay of the massive interactions, whereas at small impact parameter the 
harmonic function $(r) diverges, and as the time shift Av is proportional to $, it diverges 
as well. Other massless particles can be introduced easily: for example Maxwell interactions 
are taken into account just by having a shift (such that at Gs = we recover Rutherford 
scattering of QED): 

Gs -> Gs + g {1) g (2) /4vr (2.24) 

For transplanckian scattering, Gs ^> 1, one should take both particles as creating shock- 
waves, and these Shockwaves should interact and create a black hole. 

The generalization to higher dimensions is now pretty straightforward. Let's first notice, 
as Amati and Klimcik did also EH, that a shockwave metric 



ds 2 = -dudv + <&(x)8{u)du 2 + dx 2 (2.25) 

would shift the geodesies at u=0 by Av = $ and the S matrix was described by 't Hooft by 
the Fourier transform of the shifted wavefunction, giving essentially 

S = e iPvAv = e ip -* (2.26) 
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What we mean is that we can perform an impact parameter transform as in D=4 and get 

iA= S (S^ e ^ (el5(M) " l} (2 - 27) 

with (/i = y/s/2 = p = photon energy and p_ = /i also) 

%s)=J) <)< t (6) = _; m = aD3iD _^ (2.28) 
so a = 47t/(Qe)_3(D — 4)). Then one obtains (with q 2 = t) 

D-4 „ 

D-4 



(2n) D - 2 q D 



M = oAn- 2 n- 2 / d^(e fa '* -l)J^(z) 



.4 



dzz^- ( e ia/z - l)Jn^ (z) (2.29) 



q D ~ 2 L 



where the q dependence of the integral comes from a = aGsq "* = aGst^^ and z = qb, 
and the exponential is e lS in general, so for small 5 the bracket in the integral is i5. The 
integral can also be rewritten as 

dll 3D-8 1 

U -Wp=$)(e iau - 1)J d -a(u"b=i) (2.30) 



A-D 



but we can find no analytic expression for it. At most one can make an expansion in a which 
gives for the integral = iac, c = 2^~ D ^ 2 /T((D — 4)/2), and so 



t v (2ir) D - 2 ' nt (2tt) 



But this is an expansion in Gst^ 0-4 ^ 2 and so in D=10 we have Gst 3 <C 1, or g s (a's)(a't) 3 <C 
1, certainly satisfied. Note also that this result matches in D=4 what one obtains by ex- 
panding in Gs. 



3 Black hole production via Aichelburg-Sexl wave scat- 
tering 

Let us now analyze black hole production in the high energy collision of particles (Gs ^> 1). 
We will analyze the collision of two massless particles in flat space, in D=4 and D > 4, first 
reviewing the treatment of Eardley and Giddings ^T]. As noted by 't Hooft and analyzed 
by in this regime we have to take into account the gravitational field created by both 
particles, so one has to analyze the scattering of two A-S waves. For an estimate of the 
gravitational energy being radiated away in the high energy collision see |2E] • 
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As one can imagine, in general, the collision of two gravitational waves is a highly nonlin- 
ear and nontrivial process, and as such it is hard to say anything about the collision region. 
If we denote by I the region u < 0, v < before the collision, by II the region u > 0, v < 
(after the wave at u=0 has passed), III for u < 0, v > (after the wave at v=0 has passed) 
and IV for u > 0, v > (the interacting region, after both waves have passed), the solution in 
region IV was calculated in [21] only perturbatively in the distance away from the interaction 
point u — v — 0. 

In the case of sourceless waves (pure gravitational waves), Khan and Penrose [2H] and 
Szekeres j3H| have found complete interacting solutions, but they don't represent the collision 
of photons. We will discuss them in a next subsection. A general treatment of collision of 
pure gravitational waves can be found in j3T] , as well as in pl2"l 133] 



3.1 Review 

Coming back to the case of the collision of two A-S waves, there is an observation, first due 
to Penrose, and extended by Eardley and Giddings, which permits one to say that there will 
be a black hole in the future of the collision without actually calculating the gravitational 
field. One can prove the existence of a trapped surface, and then one knows that the future 
of the solution will involve a black hole whose horizon will be outside the trapped surface. 

An aparent horizon is the outermost marginally trapped surface. The existence of a 
marginally trapped surface thus implies an aparent horizon outside it. A marginally trapped 
surface is defined as a closed spacelike D-2 surface, the outer null normals (in both future- 
directed directions) of which have zero convergence. In physical terms, what this means 
is that there is a closed surface whose normal null geodesies (light rays) don't diverge, so 
are trapped by gravity. For a Schwarzschild black hole, the marginally trapped surface is a 
sphere around the singularity, that happens to coincide with the horizon. 

Convergence is easier to define in the case of a congruence of timelike geodesies. For a 
congruence of timelike geodesies characterized by the tangent vector £ a ,£ a £ a = — 1, defining 
B a b = Vfc£ a and the projector onto the subspace orthogonal to £ a , h ab = g ab + £ a £ b (induced 
metric), the convergence is 9 = B ab h a b. 

But we need the case of null geodesies, that is more involved. We have to first define the 
affine parameter A along the curve C such that 

Then we define a ( "pseudoorthonormal" ) basis for the tangent space, E±, E 2 , E%, such 
that El = £ a , and E$ = L a is another null vector: E 3 ■ E 3 = E 4 ■ E 4 = 0, E x ■ E x = E 2 ■ E 2 = 1 
and E± s 2 orthogonal to £3,4, but E^ b g a b = —1. And if m,n takes the values 1,2 in the above 
basis, then 9 = £, m ;ng mn - If the geodesies are null, one cannot find an orthonormal basis (as 
in the timelike case), one can only find this pseudo-orthonormal basis. Also note that by 
definition, the null geodesies defined by £ a are normal to the 2d surface spanned by Ei, E 2 , 
and we are taking the derivative of £ just in those directions. 

So to calculate the existence of a marginally trapped surface, we first need to find the 
null geodesies normal to the surface, and then impose that their convergence is zero. 
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To calculate the convergence, take the approach from j^!]. The convergence is 



6 = h ab D a £ h - £ = ?d a = —d a (3.2) 

dX 

for a congruence of null geodesies £ M normal to the surface B, and h a b is the induced metric 
on B. B is spanned by the Ei,E 2 of before, and contracting with the induced metric is 
equivalent to contracting with g mn in the above basis. 

We see that we need to get the form of the geodesies x a (X) to proceed. We can impose 
the fact that the geodesies are null, so (£, £) = 0, normal to the generators of the surface, K i: 
so (£, Ki) = 0, and also the normalization (£, dt) = —E (which can be chosen to be = —1 for 
simplicity). Note that in [31] B is a sphere, so the generators are d^. Then one calculates 
x a (X) and then £(A(x a )) = £(x a ), and then h mn = d m X a d n X b g a b (where X a are coordinates 
on the surface B), and finally 9 = h mn D m ^ n . 

Let's apply this procedure to the metric of two colliding general Shockwaves (without 
specifying for the moment the Aichelburg-Sexl solution for one moving in the u direction, 
and the other in the v direction. 

ds 2 = -dudv + dx 12 + $ 1 (^)<5(n)rfn 2 + ® 2 (x)5{v)dv 2 (3.3) 

After the coordinate transformation 

(V$ 2 ) 2 



4 

(V*i] 



it becomes 
where 



u = u + <MCw) + v9(v) 
v = v + <&i9(u) + u9(u / 

11 1) 

x l = x l + -d i ^ 1 (x)6(u) + -d^ 2 (x)9(v) (3.4) 

ds i = _ dudv + [H$H$ + H$H® - S„ d.r'd.r> (3.5) 



H\f = 5 l] + l -d l d^u9{u) 

H\f = 5ij + l -d l d 3 & 2 \9{v) (3.6) 



At zero impact parameter (b=0), and for A-S Shockwaves in D=4, we have 



$! = $ 2 = SGfilnp] p = y¥¥ (3.7) 

In general D, but for an A-S wave at b=0, there is a D-2 dimensional trapped surface 
consisting of two disks (balls), parametrized by x, of radius p c in p. 

In complete generality, the surface S is defined as follows. Take the union of the two null 
hypersurfaces v < 0, u = and u < 0, v — with a D-2 dimensional intersection u = v = 0, 
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that intersects on its turn S on a D-3 surface C. (a priori, two D-2 surfaces intersect on a 
D-4 surface though, more on that later). Then S is composed of 
"disk" 1- {v = -^i(x), u = 0}, (tf i = on C) 

"disk" 2- {u = —^2{x), v = 0} (\&2 = on C). As we will show, the condition of zero 
convergence implies that 

V 2 (^i-$i) = (3.8) 

interior to C. We will see that in the b=0 A-S case, we can actually choose = $i, ^2 — ^2 
which, with the definition 6(0) = 1, means that both disks correspond to u = v = 0. So we 
wouldn't see the topology in the bar coordinates, we need to go to the unbarred ones to get 
explicit formulas. 

On the first disk, we have 

1 u 2 

d* = -dudv + & + -u0(u) + -e^s^sfim^ 

= —dudv + dx l dx^ gij (3.9) 
and the null geodesies through {v = — u = 0} are defined by 

£ = u-?- + v£-+tfd i (3.10) 
ou ov 

The tangent generators of the surface are 

dX» 
dxi 



K j=^nr ( 3 - n ) 



where x l are the coordinates on S and the coordinates on the space, but we choose 
x % = X 1 and so 

K$ = (0, , 5)) Ktdp = -djVd v + dj (3.12) 

And we have to impose the condition that £ is null (£, £) = 0, transverse to all the generators: 
(£, Kj) = 0, and we have to define the time direction (in [31], that was (£, d t ) = —E, where 
E can be scaled to 1), in this case (£, d v ) — — 1. 
These conditions together fix 

u = 2; x i = -g ij djV; v = -d^d^g^ (3.13) 

and then we calculate 

£ = -dv - -g ij d^dj^du - d^dx 1 (3.14) 

and thus 

#| u=0 = -V 2 (^-$) (3.15) 

as advertised. 

Actually, what we have found is that by imposing (£1, d v ) = —1, we get 

& = -dv - -iy^xfdu-di^xdx 1 (3.16) 
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but similarly if we impose instead d u ) = — 1, we get 



^ = - du -jvkr- d "- i <^ dx '- (3 - 17) 

Then on disk 2, (£ 2 ,<9 U ) = — 1 implies 

£ 2 = -du-^{V<5> 2 ) 2 dv -d i V 2 dx i (3.18) 

And these two surfaces intersect on C, thus the normal, £, has to be continous across C. This 
means that for the A-S wave at b=0, when $! = $ 2 implying ^ = \I/ 2 , we need to have 

(Wi) 2 = (W 2 ) 2 = 4 (3.19) 

Then in D=4, replacing the explicit form of <j) we get 

$ = $ = SGfi In p/p c p c = AGfi = r h (3.20) 

whereas for D > 4 

<3 - 2i) 

In the bar coordinates, both disks correspond as we said to u — v = and x % = x l . But 
this surface in the bar coordinates is just fiat (on it, the metric is Minkowski), so the area 
(volume of balls) is just the area of two flat balls of radius p c = p c . The area (volume) of a 
flat unit D dimensional ball is Ha«i,D = &>d-i/D, so the total area of the trapped surface in 
D spacetime dimensions (two flat balls) is 

A mm (S) = 2V ball , D - 2 p?- 2 = £^^-3P?- 2 (3.22) 

whereas, from the explicit form of the Schwarzschild solution in D dimensions the horizon 
radius of a black hole of mass y/s = 2/x is 

so that the horizon of the mass= \fs black hole is 

A SA = IWf" =► = i[ (0 - 2)n " ]A ^ \ (3.24) 

The area of the trapped surface is smaller than the horizon area of the black hole to form 
(since the horizon is by definition outside the trapped surface), and we can express the area 
of the disks (balls) as the area of horizon spheres that will form, so r < r^, where r is defined 
by Area(S)= Oo_ 2 r D ~ 2 , implying that the mass of the formed black hole satisfies 

IgTTgMgg n_ s > r D-3 = ^rea(S) ^3 Mm 1 (g - 2)n g . 2 _i_ 
(L>-2)fi c _ 2 h ~ ~ [ Q D _ 2 J "2 [ 2fi D _3 J l ' J 



(we have put in the explicit form of Area(S) and of p c in terms of \x = yfs/2). Both 
and A '^ g match the explicit numbers in 
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a. 



3.2 Extension 



In the previous discussion we have already set up the formalism so that it is valid for any 
function characterizing the shockwave. We will be applying this later for different $'s. 

Let us now try to extend this for the case of nonzero b in any dimension. For the b=0, 
D=4 A-S wave we had 

■qr = $ = -8G/1 In p/p c (3.26) 
meaning that $f > for p < p c . For b=0, D > 4 we have 

* = *_* (p = „ o); * = _J^_ (3 ,T) 

and again \I/ > for p < p c . 

For b > 0, D > 4 now, we would need both ipi and ^2 to be zero on the same surface 
(curve, for D=4) C, not on two surfaces C\ and C2, since then the intersection of C\ and C2 
would have D-4 dimensions (points, for D=4). So we cannot use in D=4 for instance 

^ = = -8Gp In ' P ~ Po ^ (3.28) 

That would define two disks in x = x with the centers dispaced by b, and while each of the 
disks boundary would be a circle, the two circles will intersect in two points, so C would be 
composed of these two points. 

The correct solution, which was explored in using a self-consistent approach (which 
does not guarantee to find ALL solutions) is that ^ 7^ $ i; just 

V 2^. = V 2 $ . x _ - o .) ^ 3 2g ^ 

which means that ^1, ^2 are Green's functions for sources at xqi, X02 which both are zero on 
the same curve C enclosing xoi and xq2- Then one imposes the condition for continuity of 
the null normal £ which gives 

■ W 2 = 4 (3.30) 

which fixes (together with the previous conditions) the form of C. 

Clearly for very small b (much smaller than p c ) we have that C is well approximated by 
the boundary (envelope) of the union of the two disks. We will assume that in D=4 the two 
points C\ and C2 (intersection of the two circles=boundaries of disks) are still part of the 
curve C even at b large, which seems like a reasonable assumption, though not well justified. 
Let's see what we can deduce out of it. Clearly the two sources will be inside C, so we will 
therefore assume that the curve C is outside the paralelogram made from Ci, C 2 , x i, Xq 2 - 
We still call the distance between Ci^ and £01, £02 ( ra dius of the circles) p c , and if we then 
impose (J3.30|) on = $1 and ^2 = $2 (which are still good Green's functions for the circles 
that both pass through the two points C\,C2, even if they are not for the whole curve C) 



we get the equation (cos 8/2 = \/ p 2 c — b 2 /A/p 



vl/i \1>9 R 2 b 2 

|Vf M Vf|^ = l^-(l-^) = l (3.31) 
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(where R s = 4/iG) which gives the value of p c as 



Pl = -f^ + \j l - ( 3 - 32 ) 

We can check that if b=0 we reproduce the known result of p c = R s . This formula means 
that the maximum impact parameter for which we can have a black hole forming within this 
approximate formalism is b max = R s /\/2 = 4G/i/v / 2 and the minimum radius is p c ,min = 
Pc(bmax) = R s /V2 = b max , and the area of the trapped surface satisfies 




2b 2 b 2 

- #r) - 2" = Smm (3 - 33) 

so that S m in at the maximum b is S m in — -R^V / 3/4 = Ay/3(pG) 2 . 

Comparing now with the results of jTTj we have that b max = AGp/ \/2 ~ 2.83G/X is smaller 
than their result of 3.219GyU. Since b max < R s it is even physically acceptable (we would 
have a problem if it would be bigger). As for the estimate of the area of the trapped surface, 
Smin — 4v^3(/iG) 2 , it is sensibly smaller than the result of which can be found to be 
(replacing the value of their parameter a max in the formula for the area) 40.852(/iG) 2 , so we 
have a much more conservative estimate. 

But the advantage is that this procedure can be now easily extended in higher dimensions. 

Indeed, in D > 4 , 

\QnjiG x 

&D-3 P 1 



and so the condition V$i • V$2 — 4 implies 



2 



(^)-(l - ^ = 1 (3,5) 

where R s = is the horizon radius of the black hole, e is defined in f!3.24j) and this equation 
can be rewritten as 

f{x) = Ax D ~ 2 - Aax + 2b 2 a = (3.36) 

where x = p 2 . and a = (eR s ) 2D ~ e . We can easily find the maximum value of the impact 
parameter b from it. Since p c = \fx is the biggest of the solutions to the equation ()3.3(ij) . we 
impose that f(xo) < 0, where xq is the highest root of f'(xo) = 0. This condition implies 

,9 r a . i D-3 2(eR s ) 2 , . . 

b 2 < 2 — — u=i— — = 1 S D _ 2 {D - 3) (3.37) 
D-2 D-2 [£)-2]s=s 

We can check that indeed in D=4 we recover the result b max = R s /^/2, since then e = 1. In 
D=5, that means b < 0.9523i? s < R s . We can also calculate the lower limit on the area of 
the trapped surface as before, except that now the area of the paralelogram Ci, C2, X01, X02 
is replaced by the volume of a "surface of revolution" in D-4 transverse directions around 
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the axis xo\,xo2- The geometry in higher dimensions is more complicated, but for D=5, this 
is just two cones glued on their bases, of height h=b/2 and base radius p c cos9/2 , and so 
"Smin" (volume of the cones) is 

Snh ixb , 9 b 2 , . 
2^- = y(p^- T ) (3.38) 

In conclusion, we have set up a formalism for Shockwaves in which we can calculate 
trapped surfaces at b=0 and to some degree at nonzero b, for a general shockwave form. 



3.3 Collision of sourceless waves 

We have seen that for A-S-type waves colliding, in general we get a trapped surface in the 
future of the collision, which indicates a black hole horizon being formed. From this, we 
conclude that a black hole is formed in the high energy collision of two high energy photons 
(massless particles, with an energy momentum source). 

But what happens if two sourceless waves (gravitational solutions to the pure Einstein's 
equations) collide? We would expect to be able to associate this phenomenon with the 
collision of two gravitons, in which case we would expect to create a black hole in the 
collision. It is in fact true that there is a theorem stating that a singularity will form in the 
future of a collision of two sourceless waves |3*2*l 13*3*] . It is also a theorem that for Einstein 
gravity in flat space, a singularity cannot be naked, so we would expect to be able to find a 
trapped surface, indicating the formation of a horizon in a sourceless wave collision. 

Unfortunately, we will see that this is not so, and we will speculate on why, after we see 
the problem. 

Khan and Penrose found a solution j2H] describing the (head-on, zero impact parameter) 
collision of two source-free gravitational pp waves of the type 

ds 2 = -dU(2dV + (X 2 - Y 2 )h{U)dU) + dX 2 + dY 2 (3.39) 

with h(U) = S(U) (the function <3> defined before satisfies the source-free equation <9,f$ = 
solved by $ = —X 2 + Y 2 ). After the coordinate transformation 

U = u, V = v + x 2 /2FF' + y 2 /2GG'; X = xF; Y = yG (3.40) 

withF" = —Fh, G" = Gh, solved by F = l-u9(u), G = l+u9{u) (F' = 9(u), as u8{u) = 0, 
and thus also F5(u) = S(u)), we get the wave in the form 

ds 2 = -2dudv + F 2 dx 2 + G 2 dy 2 (3.41) 

The collision will involve two such waves, one in u and the other in v, at zero impact 
parameter (b). Thus the colliding wave solution of Khan and Penrose is 

ds = — ^ + t{ )( )dx +t ( — — — ) ( — — — (3.42) 

rw{pq + rw) z r — q w — p r + q w+p 
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where 



p = u9(u); q = v9(v); r = a/1 — p 2 \ w = a/1 — q 2 ; t = a/1 — p 2 — q 2 (3.43) 

In the region u > 0, v < (before the coming of the second wave), we can check that the 
Khan-Penrose solution becomes 

ds 2 = -2dudv + (1 + pfdx 2 + (1 - pfdy 2 (3.44) 



that is, of the sourceless wave form (j3.41j) . and we see that there is a coordinate singularity 
at u=l. Then in the collision region u > 0, v > we have 




(y/(l-U*)(l-V*)-Uv) 

dudv 



){^{l-u 2 ){l-v 2 ) + uv) 



u 



\ 2 dx 2 + (VT^u 2 - vfiVT^v 2 - ufdy 2 



u 2 — V 2 



(3.45) 



Putting v=0 we get back to the sourceless wave solution ()3.41j) . And [2H] found that in the 
collision region, the line u 2 + v 2 = 1 has a scalar curvature singularity. We can calculate that 
for u 2 + v 2 = 1 — e the metric is 



yfedu 2 Auv s 
2u 2 v 4 + ^~T' "™ ' y 4uv' 



^ = e[-H-r + ( — ) dx + br-rdy ] ( 3 - 46 ) 



so clearly the metric is singular, but there doesn't seem to be any good way to define a finite 
area of the singularity. Indeed, at u=fixed, v=fixed, 

dS = ds x ■ ds y = edxdy (3.47) 

So we can't calculate this way a minimum on a horizon area of a black hole that would 
probably form. 

But we can still try to apply the formalism of Eardley and Giddings and calculate the area 
of a trapped surface that we assume will (should) form. Indeed, the individual gravitational 
waves that collide are still of the general form used in the previous subsection. The only 
difference is that instead of $ = — 8Gp, In p/p c we have 

$ = - X2 ~ y2 ; (p c = 4G/i) (3.48) 

Pc 

where we have rescaled U and V to introduce the dimensionful parameter p c describing the 
strength of the wave. The problem is though that in order to be able to choose \P = $ 
like we did for the b=0 A-S wave collision, $ would have to be zero on the curve C at its 
boundary, so that p = p c as we shall see. But $ = for X = ±Y (we could shift $ by a 
constant, and then C would be a hyperbola), so that we can't actually choose ^ = $. 

Thus we could only use the above function for X = ±Y, and these are four points 
(X = ±p c / y/2, X = ±Y) that would presumably lie on the curve C if there would be such 
a curve. 
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However, the correct treatment would involve solving the 2d Green's function ( "electric 
potential") for the Laplace equation V 2 \l/ = with Dirichlet boundary conditions \1/ = on 
a curve C where V\& (the "electric field") has unit norm, and this is impossible! 

Thus we seem to have proven that the assumption of a trapped surface is in fact wrong! 

So there really seems to be no way of obtaining a trapped surface in the Khan-Penrose 
solution, even though we do obtain a singularity! 

This is a most bizzare situation in itself, which could be perhaps saved by the fact that 
in such singular spacetimes the usual censorship theorems don't apply, but correlated with 
the expectation that the Khan-Penrose solution should describe graviton-graviton scattering, 
this is really puzzling. 

One could perhaps think that the Khan-Penrose metric is not the correct sourceless wave 
to describe graviton-graviton scattering. After all, there is a plethora of sourceless wave 
scattering solutions, as reviewed to certain extent in [3T] . 

One can analyze their behaviour though (we will not do it here explicitly) and convince 
oneself that these solutions do not describe graviton scattering. The simplest of them is 
the Szekeres solution in jSU], which has the same singularity structure as the Khan-Penrose 
solution, but is described by a function $ of the type $(m) = f(u)6(u) as opposed to the 
delta function profile $(£t) = 5(u) of the incoming waves in the Khan-Penrose solution. 
The rest of the possible solutions are even more complicated, and they really describe the 
collision of realistic gravitational waves, as opposed to the idealized delta function waves of 
the Khan-Penrose solution. Therefore the Khan-Penrose solution is the only one that can 
claim to represent the collision of two (idealized) gravitons. 

Of course, all these solutions were in 4d general relativity. 

Gutperle and Pioline [23] set out to generalize these solutions to 2n+2 dimensions and to 
add p-form field strength to it, the ultimate goal being to scatter 2 maximally susy pp waves 
of IIB, or rather a shockwave generalization of it. They fall kind of short of the goal. The 
first try at the generalization gives exact solutions which however do not satisfy appropriate 
boundary conditions: the incoming waves are different from the Khan-Penrose and Szekeres 
profiles. 

A perturbative attempt near the lightcone (or for the strength of a wave much smaller 
than the other) produces a higher dimensional solution, as well as the p-form generalization. 

Then, j2H] also produce some generalizations of this type (see also [23 EH]), with better 
singularity structure, but they don't analyze the incoming waves in Brinkman form, so it is 
not clear what they correspond to. 

In conclusion, the Khan-Penrose solution is the only one that has a chance of describing 
the collision of two idealized gravitons, and we seem to obtain the existence of a naked 
singularity (no black hole!) in the future of the collision. A good explanation of this paradox 
is still lacking. 

4 String corrections 

We will now try to apply the previously derived formalism to shockwave metrics that incor- 
porate string corrections to the high-energy scattering of two photons. 
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There are two such formalisms. One is due to Amati and Klimcik [21] , and the other due 
to Amati, Ciafaloni and Veneziano [22] (see also [3§] HO] l4*T]). 

The approach by [2H] involves writing down an effective shockwave metric from which one 
can calculate an S matrix, which then is matched with a string-corrected S matrix (string 
calculation). 

The S matrix is defined as exp(iS/h) where the action is a function of the classical 
effective metric, with a source coupling to an external T^ u . Namely, the S matrix is 

S eff (b, E) =< e iA{h ^ )/h > tree = e tA[h ^ )/h (4.1) 

where 

A(hT) = J d 4 x[L eff (h^) + T, W W V \ (4.2) 

is the action evaluated on its classical solution with sources given by two Shockwaves at x=0 
and x=b 

T__ = kE5(x~)5 2 (x); T ++ = kES(x + )5 2 (x - b) (4.3) 

and L e jf is an effective Lagrangian by Lipatov |42| . Amati et al. [2S] showed as we mentioned 
that this calculation reproduces the result for the string correction to the scattering matrix 
S. 

The string-corrected A-S type metric obtained in j2H] for D=4 can be expressed in terms 
of a $ of the form 

$ = - S) + = kE[-—log^r + R 2 s a {l) (z)} (4.4) 

where 

and kE = SnGfi (k = SiiG,E = /i), so the coefficient of the log term in <3>(°) is —R s . Here 
z = x 1 + ix 2 are complex transverse coordinates so that \b — z\ 2 — (b — X\) 2 + x 2 , etc. 
Then 

= h{x *l xl) M -A- bx 1 )log {b - X ^ + X2 + X f(x 2 + x 2 2 - 2te0] 

= u - b)log {h - X ^ + xl + \{x\ + x\ - 26x0] (4-6) 

b[x\ + X2) b 

But the problem is that the string-corrected metric is only valid for b > R s (when no black 
hole forms yet), whereas we want to have a perturbation in b small. 

We have tried to just plug in this metric in the continuity condition V^i ■ V^2 = 4, and 
treat it perturbatively in R s /b as in j2H], but one gets corrections of the order R 2 /p 2 c1 and 
then in a perturbative solution, when p c is replaced by the first order solution which is of 
o(R s ), the corrections are of order one. Thus this perturbation is useless. 
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But one can still do a small calculation, namely to take the corrected metric (with b 
nonzero and moreover > R s ) and see what it does to the continuity condition for the collision 
at b=0 of two non-corrected metrics, namely (V$) 2 = 4, now with 

d^ = -2R s ^- 2 +dM 1] (4.7) 
P 

and expand in R 2 s /b 2 . This will not be very relevant (since the corrections dissappear in 
the b — > oo limit, which is not what we want), but just to see what kind of effect string 
corrections have. We obtain 

p\ = R\[l -^A + - xl)lo S Y l ~ ^' + - ) + A] 

A = + („ - b )lo 9 &^±^ ) (4.8) 

b er 

where we should have x\ + x\ = p 2 , but obviously since we have used an asymmetric solution 
(where b is a distance on the X\ axis), the solution we get for p c (R s ) also depends of our 
choice for x±,X2, so p c = p c (xi,X2). 

The above solution is exact, but we only need to expand in p c jb. Expanding to the first 
two nontrivial orders we get (after some algebra) 

p\ * ifll + + x\) + ±(8xlxl -Rj + Ax\R 2 s - \x\)..] (4.9) 

The first correction is proportional to X\ (times a positive quantity), so when we calculate 
the area of the curve p c (xi,X2), the positive contribution for x\ > will cancel against the 
negative one for x% < 0. So we need to turn to the next correction to see whether or not the 
area increases. 
Defining 

f(xl) = Sxlxl -R 4 S + 4xlR 2 s - % -x\ ~ \2x\R 2 s - ^x{ - R 4 S (4.10) 

for y = xf/R 2 between and 1, we can check that the function is positive for y > 0.09 (most 
of the domain), so a simple estimate shows that the area of the trapped surface will indeed 
increase. 

But after so many approximations it is not clear this still is relevant. 
We turn instead to the approach of Amati and Klimcik [21] . 

Amati and Klimcik [21] first generalize the 't Hooft and Dray and 't Hooft calculation, 
as we explained in section 2. A shockwave metric 

ds 2 = -dudv + $(x)5(u)du 2 + dx 2 (4.11) 

would shift the geodesies at u=0 by Av = $ and the S matrix was described by 't Hooft by 
the Fourier transform of the shifted wavefunction, giving essentially 

S = e ip " Av (4.12) 
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In string theory, the 't Hooft scattering in the shockwave background gives (for an open 
string—* photon) 

Av= - I ${X{(T,0))da (4.13) 

Jo 

and the S matrix is defined as acting on creation/annihilation operators as S + cii n S = a ou t. 
Then 

This matches the resummed string calculation of jlU] if 



*(!/) = "ff" r - ■ *tree(s, y - X d (a d , 0)) : ^ (4.15) 

Jo S IT 

where 2pvq = — 1, b = x u — x d and X u , X d are nonzero-modes. Here the indices u,d refer to 

s ' ' ' 

"up" and "down", necessary when we evaluate $(X(cr, 0)). 

We note that here b refers just to a parameter in the calculation of the shape of one 
modified A-S metric. We haven't reached the scattering of two A-S type waves yet. In that 
case, we will denote the impact parameter of the two waves by B, to avoid confusion. 

Then we match with the S matrix obtained by resumming string diagrams, 

S = exp[2i / : a tree (s, b + X u (a u , 0) - X(a d , 0)) : — V^] (4.16) 
Jo 7r 

and the tree amplitude is 

a ^ee(s, b) = ^z- 2 b'~ D / dte-H D '^ (4.17) 



(Y = a! logs). Then $(?/) becomes the function for the A-S wave at y ^> yY = y/ a' logs. 
So S is dominated by graviton exchange at large b (Aichelburg-Sexl) and by absorbtion at 
small b. 

As a first approximation, we can neglect all string oscillators in $(y) and obtain 

= a tree (s,x) (4.18) 

s 

where a tree {s, x) is a tre e{s, b) (impact parameter space), and becomes equal to the A-S result 
at large b. We can rewrite it also as (g= gauge coupling) 

a, (b s) = - C dpp D / 2 " 3 e~wP = t^L I ^ dte-H D / 2 - 3 

(4.19) 

where Y = a'log(-is) =Y — ina'/2. We are only interested in the real part of a tree , as it is 
the only one that we can use in the classical gravitational wave scattering calculation. It is 
obtained jut by replacing Y with Y. For b 2 3> Y, we obtain 

Rea tree (s,b) - f ^ D/2 ^4 (rW2 ~ 2) - e -^(^)^ (1 + {D/2 _ 3)^) + ..,) 
q 2 sa' , 2 . . 
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whereas for & 2 C F we get 

g 2 sa' 2 1 b 2 

Re a tree (s, b) * ^J^yy^i^ - AY(D _ 2) + -) (4-21) 

One need just repeat the Eardley-Giddings-type calculation now, as we have set it up in the 
previous section. 

The regime we are working in is small g, large G^s = g 2 a's/(8ir). Since R 2 = AG 2 N s = 
g A a' 2 s/(A-K) 2 and Y = a'log(a's), 

Rs 9 2 g 2 ois 

~Y = log{a's) (4vr) 2 ^ ' ' 

can still be arbitrary, in particular it can be very large. Since to first order b max = p c = R s (for 
Aichelburg-Sexl), and at large b the metric is A-S plus corrections, in the regime R 2 /Y 3> 1 
we can use the large b (b 2 /Y ^> 1) expansion of &(b). 

Then in D=4 we get, with p u q v = s =^> q v = \/s (with the choice p u = q v due to center of 
mass scattering, with equal strength Shockwaves scattering) 

*(b) = - 9 ^a> (2 log±--e-^(^)- 1 + ...) = -R s (2 log ±- - e"^)" 1 + ■■■) ( 4 - 2 3) 

Then the condition for the trapped surface appearing in the scattering of two Amati-Klimcik 
waves at zero impact parameter, (V$) 2 = 4 gives 

AY Rj 

b max = Pc^ R s (l + (1 + ~5?) e ~ ) (4-24) 

Rj 

(for b 2 /Y ^> 1, so R 2 /Y ^> 1) thus increases, so the area of the formed black hole also 
increases (since the black hole area is proportional to p 2 ) . The area of the trapped surface 
giving the bound on the horizon area is SVmn = = an d r h = 2MbhG, so 

M " = 57fe (4 ' 25) 

also increases. 

At nonzero impact parameter of the two Amati-Klimcik waves, parameter denoted by B 
as we mentioned (to avoid confusion with the b that was used previously), applying the same 
approximation for finding p c as was used in the flat space A-S case, the normal continuity 
condition is ■ d$2 — 4, so <9j$(x — X\) ■ <9j$(x — x 2 ) = 4, so we only get an extra factor 
of 

B 2 



cos 2 6 = 1 -— (4.26) 
to the condition, which thus gets modified to 



2PI 



f=v i -5 (i+e_4+ - ) (4 ' 27) 
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solved perturbatively by 




_1„ , /, 2 B2|o; „ + yo= ^Vlll (4.28) 



which means that £ maa; = R s /y/2{\ + e- fl ' /(8y) ). 

Finally, let us see what happens if R 2 S /Y <C 1. At first, we would guess that we can use 
the small b expansion of the metric b 2 /Y <C 1, for which 

But if we plug it into the continuity equation for getting p c , (V$) 2 = 4, we would get 
p c = 4Y/R S to first order, meaning that p 2 /(4Y) = 4Y/R 2 , that is we would seem to be 
in the opposite regime, so the perturbation expansion used was invalid! The solution is 
of course that R 2 S /Y <C 1 will correspond to p 2 /Y ~ 1, so we would need to use the full 
solution, which however is difficult to handle. 

But in any case we can say that for R 2 S /Y <C 1, classically (A-S wave) we have p c = R s , 
but in string theory we get p c ~ y/Y ^> R s , so we have a great increase in the area of the 



black hole formed, thus it is natural to assume the cross section will also increase. 



5 Randall- Sundrum- type models 

The next application of the black hole creation formalism is to see what kind of correc- 
tions appear if we have the black hole being created in a physical setting, namely for a 
Randall-Sundrum scenario for low Planck scale. Emparan j33] found an A-S-type wave in 
the background of the one brane RS scenario, and analyzed the scattering a la 't Hooft in 
this wave. 

Here we will try to see how the addition of the RS background affects the Eardley- 
Giddings calculation for the flat space black hole creation. We will keep the wave on the 
brane, as in the Emparan calculation. 

5.1 A first attempt- applying the formalism 

The solution for an A-S-type wave in the RS background is 

ds 2 = dy 2 + e~ 2lvl/l (-dudv + dx i dx i + h uu (u, x\ y)du 2 ) (5.1) 

where 

Ku = ( 27r )(M/2 ^( M W J dqq VWgr) Kd/2 _ i{lq) (5-2) 
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which is a solution of Einstein's equation with t uu = 2irp5{qo + q±). Yet another form for the 
metric is 



e- 2 ^ l h uu (u,r,y) = -4G,p6{u)[e- 2 ^ l log- 



21 [°° .Y 1 {ml)J 2 {mle^ 1 ) - J^ml^mle^/ 1 ). 
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dmK ° {mr) J 2 {ml)+Y 2 {ml) ] (5 ' 3) 



which means that on the brane (y=0) 

h uu (u,r,y = 0) = -AG 4 pS{u)[log^ - A f° — ] (5.4) 

I 2 n Jo m J ({ml) + Y^{ml) 

The Einstein tensor for this solution is linear in h uu , and thus even though this is found as 
a solution to the linearized equations of motion, it is also an exact solution. 
At large distances, r 3> I, 

r 2 I 2 2/ 4 r 2 
h uu {u, r;y = 0) = -AG A p5{u)[log— - — + —(log^ - 1) + ...] (5.5) 

whereas at small distances r <C I, 

I 3 v 3?^ 

h uu {u,r;y = 0) = -AG A p5{u)[-- + -log- + — + ...] (5.6) 

We can use the formalism developped previously, since the solution can also be expressed 
as just a modification of the $ function. Now we can at least calculate the zero impact pa- 
rameter (b) values of S min (the area of the trapped surface) and the mass of the corresponding 
black hole. We can also estimate the nonzero b parameter values of p c {R s ), b max , S min . 

The new function $ is now 

P 2 4 f°° dm K (mp) 

*(»,p,y = °) = -Wo4-- 2 j o - JMi ;\^ (ml) ] (5.7) 

which means that 

= - ± r dm ffi^ i (5 . 8) 

P P tt 2 Jo J ({ml) + Y{{mly 

and thus imposing the continuity of the normal condition (<9j<3>) 2 = 4 and rescaling the 
variables by R s we get the integral equation for p c 

Pc _ 1 2p e /R. f 00 K{yp c /R s ) 

Rs * 2 Jo y j 2 {yl/R s )+Y 2 {yl/R s ) 

As before, the area of the trapped surface is the area of two disks, so it is 

S min = 27T P 2 c = A7rr 2 h (5.10) 
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where is the horizon radius of the formed black hole, and = 2GMj,h, so 



We can use the expansion for p ^> / and p <C I to calculate the form of p c from the 
equation (|5.9|) . for i? s 3> / and R s I. For R s ^> I we have 

$ = $(0) + $(D ; $(l)^ jRs £_^ ( ^_ l) + ... ] (5 . 12 ) 

and thus imposing (<9i$) 2 = 4 we get 

p 2 ~ R 2 S [1 + ^ - ij(8 Zn|f - 13)] (5.13) 



For i? s <S / we get 



and then 



• = -*[- + §"•? + ! + ■■■] (514) 




p c ~ W— i(l + -W_f + --f + ...) (5.15) 



Note that p c = R s is what one gets in flat 4 dimensions, whereas p c = y/2G^p = y R s l/2 
is what one gets in flat 5 dimensions, so the formula is correct to zero-th order. 
So the mass of the black hole is 

M 6 ,~^(l + |l + ...) Kfl. 

^ ^ ^\Z3r (1+ tVW + "° />>i?s (5 - 16) 

Notice that the limit of small 1 is the limit in which the space is very 4-dimensional (large 
exponential warping in the extra dimension), so the four-dimensional result should hold, and 
we find that (just small corrections to the usual 4d result). The limit of large 1 is when the 
background space is approximately flat 5d space, so we have to modify the results to account 
for the creation of a 5d black hole. The condition (V$) 2 = 4 is independent of dimension, 
but it becomes (<9j$) 2 + (d y Q) 2 = 4 in a general dimension (with y being the transverse 
dimensions), and it will be modified for a general background. 

Thus in the general case the trapped surface is something in between two disks and 2 
balls, so two fat disks, or flattened balls. In the 2 limiting cases, the trapped surface can be 
approximated by 2 disks or 2 balls, respectively. One can still define the black hole projected 
onto 4 dimensions. 

We will come back to the correct treatment in the next subsection, and we will see that 
whereas the zero-th order formulas are correct, the first order corrections get modified. 
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At nonzero b, applying the same approximation for finding p c as was used in the flat 
case, the normal continuity condition <9j$i ■ <9j$2 = 4 becomes d^(x — X\) ■ d^(x — x-z) = 4, 
so we only get an extra factor of 

We 

to the condition, so that now 



cos# = l-— (5.17) 



R s V 2 P* 11 J 2 a (yi/Rs) + Y?(yi/R s y { - ] 

whereas the expression for the (very conservative) estimate of the trapped area, S m i n , remains 
the same as a function of p c and b, 



b 2 P 2 c - j (5.19) 



Expanding in the I -C R s regime we get 



o 2 b 2 2l 2 I 4 o 2 

%r {1 -M )ll+ *-7J m f- n)+ -- ] (5 - 20) 



so that 



pi 1, \ 2b 2 8l 2 , b 2 N J\, , 

|=2< 1 + V 1 -fi; + W (1 -^ ) + < ' ( ^ )) (5 ' 21) 



(so that bl ax ~ J RJ/2(1 + 4/7.R2)). 
In the / 3> R s regime we have 



IRs b 2 3 ft 3^ 

: Tf"^ 1 + 2l + 8P 



P C 2 = ^A 
The first term gives the equation 



X =a(x-j); x = p c ; a = — ( 5 - 23 ) 

Solving this equation and selecting the solution that gives x = a in the limit of b=0, we get 
(also calculating the first two corrections) 



M = ^ = a( i+ 3 + + (5.24) 

IR S a 2 V 2/ 8 2/ ; v ; 

where 

a = ^(A + -); A = { -p + V=iTF)* ; P = (5.25) 
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If f3 < 1, then a is real and is 



cos^/3 where cog ^ = A = e m (526) 

a/3/2 M V 7 



If /3 > 1, the solution is complex, thus 

«L. = ^ (5.2T) 

5.2 Correct treatment: generalizing the formalism 
to curved higher dimensional background 

Let us try to understand what happens to the black hole area when we have a curved 
spacetime background of the RS type: 

ds 2 = e~ 2lvl/l [-dudv + dx 2 } + dy 2 (5.28) 

Let us denote e~ 2 ' y " 1 = A and g^ = Ag^ represents the metric in both x and y coordinates 
(transverse). Then a straightforward calculation along the lines of the flat space case finds 
the vector normal to the surface is 

£ = --tfidi^dj^du -dv- d i ^dx i (5.29) 

and so similarly to the flat case the continuity condition for the normal is 

g ij dS>d^ = 4 (W) 2 + A(dy^) 2 = 4 (5.30) 

(the relation fixing the boundary of the trapped surface, or its radius) 

For these and the next relations it is necessary to calculate the coordinate transformation 
from the coordinate system 

ds 2 = e~ mll {-dudv + dx 2 + h5(u)du 2 ) + dy 2 (5.31) 

to the coordinate system without delta function discontinuities, up to order u (near u=0). 
The calculation is a straightforward but tedious generalization of the flat space case, and 
one finds after the coordinate transformation 

u = u 

v = v + h0{u) + —^-{dihdjhg^ + A{d y h) 2 ) 

■ , u8(u) , 

x — x -\ —g^djh 

2 



V = V + ^Ad y h (5.32) 
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that 



ds 2 = A[— dudv + dx 2 + u6{u)didjhdx l dx^] 

+dy 2 [l + uO(u)Ad 2 y h} + dydx l u6(u)Adid y h + dydAu9(u)d y h + o(u 2 ) (5.33) 

where 

A = e~ m/l + o(u) 2 A\ u=0 = e- 2Wl ; dA\ u=0 = -jA[dy + ^d y hdu] (5.34) 

The convergence of the normals 9 = g^D^j is now again 

9 = -V 2 (* -h) (5.35) 

where h uu = hS(u)(= $6(u)) and 



V 2 = -V 2 x + d 2 - lS gn(y)d y (5.36) 



iy 2 + d 2 - - 
Therefore we write 

$ = $ + V 2 C = (5.37) 

So now the trapped surface is a surface f{p,y) = defined by both = C (const) and 
by g^di^dj 1 ® = 4. In the flat case the first implied p = Po and the second po — Rs- But we 
also saw that the nonzero b case had the same problem as we have now: find a surface C 
and a function ( that satisfies both \1/ = const, and V 2 \I/ = 4 with \1/ = $ + £. 

In general it is a hard problem, but at least perturbatively, in the two limits I — > and 
I — ► oo we expect to find approximate disks and approximate balls, respectively (and fat 
disks in between). We would also expect that in the / — > the surface is the same disk 
p = R s as for flat 4d space. 

The formula for $ (h) at nonzero y is (in j^H], it's not $ but $e -2 ' y '/ z ), so 

<$> = -R s [log- - -e 2 ^ 1 dmK (mr) U j 21 ^ / 21 >-] 5.38 

/ 2 7r J J({ml) + Y{{ml) 

(and actually, this is defined up to a constant, so the log r 2 /l 2 is conventional, we could have 
log r 2 /r 2 ). 

Then we have 

8 , 2R 4 f°° d(ml) K (mr) 

y ly=0 7i [ hj ml J 2 (ml)+Y 2 {ml) 

2 [°° „ lW , , ,Y 1 (ml)J 2 (ml) - JAml)Y 2 (ml) , 

/ Jo Jf(mi) + Yf (mZ) 

where we have used that Y u (x)J' u+1 (x) — J v Yl +1 (x) = —1{y + l)/7rx 2 , which we can easily 
deduce from the Bessel function properties. 
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Then we find 

2R S 8 f°° d{ml) K (mr 



d y^ y ~° tt2 [ p J Q m i J*(rrd)+Y?(ml) 
+ H 2 y o d(ml)K (mr) ^ - ^ ( J } ^ (5.40) 

Let us now analyze the perturbation m l/r (the space is approximately flat 4d) 
Using the relation 

Y v {x)J'l +1 {x) - J,(x)F; +1 (x) = —(- 2 - 1) (5.41) 

7TX X Z 

which can be easily derived, and also the expansion 

Ji(x)~x/2; nYAx) ~ + xlog- + ... (5.42) 

x 2 



we find 

ARJ 2 



r 

We also have 



d 2 y $\ y =o = -—^ + o(l 4 y) (5.43) 



Let us expand ( near y=0 as 

,2 



r Z 2 

$| y=0 = -2RJog- + R s - + ... (5.44) 



C = Co(r) + Ci(r)2/+yC2(r) (5.45) 

Then at y=0 V 2 C = implies 

3Co(r)+C 2 (r)-yCi(r) = (5.46) 

and we don't want to upset the flat space solution, so we will take ( = (otherwise the 
continuity condition (V\I/) 2 = 4 implies a different radius for the trapped disks). So (2 = jCi- 
From 

(a^) 2 + e~ 2|? ' l// (^*) 2 =4 (5.47) 

we see that if d y ^f has a y-independent piece, we will change the continuity equation at y=0, 
and we don't want that to happen to leading order in 1. As <9$ | y=0 = already, we must put 
Ci = to leading order, so at least Ci ~ o(l), which implies (2 = o(l) as well. 
Then 

# = / + ay+y<7+... (5.48) 



where 



/ = $^=0 = -2R s logr/l + i?,/ 2 /r 2 + a = Ci 

= flS*|y=o + jCi = + ^70 + jCi ^0 + ~G + ... (5.49) 
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We have to check now that the two surfaces in (r,y) defined by \1/ = const, and the normal 
continuity are the same to first nontrivial order in y and 1. 

V 2 

= C = f + ay + ^-g+ ... (5.50) 



2 



and the other 



C = A=(f + yd + |V + ...) 2 + (1 - f + 2^ + ...)(a + yg + ...f 



= f 2 + a 2 + y(2a'f - 2j + lag) + ... (5.51) 
if a is nonzero and 

C> = 4 = f 2 + y 2 (fg> + g 2 ) + ... (5.52) 
if a=0. If a=0, we get to order y 2 (first nontrivial) for \& = C 

r I 2 2R I 2 

2RJog- - R s - 2 + ... = y 2 (-^-) (5.53) 

(we have traded C for tq) and for the continuity equation 

(4 - f 2 - a 2 =)4 - 4^(1 + 2^) -a 2 = y\g 2 + f'g') = y 2 {-^-) (5.54) 

Notice that at 1=0 the l.h.s. of the two equations would be 2R s 5r/ro and 86r/R s respectively, 
so with r = R s (from y=0) the two equations are not the same. So we have to put a nonzero 

Ci- 

Also note that since the constant C (and hence r ) is an arbitrary constant, at y=0 but 
1 nonzero we don't need to have the same 1 dependence in the two equations, we can absorb 
the unwanted 1 dependence in the redefinition of tq. The 1 dependence of the radius r max is 
deduced from the continuity equation (which doesn't have a free parameter). 

Also note that a priori one could check the values for $ and its y derivatives by using the 
alternative solution for $ in [3^] . We have tried to use perturbation theory on the alternate 
form (integral of ratio of K function) of $, but as Emparan noted, it is much harder to do so. 
In particular, one has to use the freedom to add an arbitrary constant to h (this is related 
to a rescaling of u and v). 

If now we put a = £i ^ (and so g = g + y 1 ), the first order in y is linear, and by 
requiring that at 1=0 we get the same y dependence in both equations we get the condition 

a^^(2ag + ^-^a>) (5.55) 
4 / r 

Thus if we put 

a = — (5.56) 
at 1=0 and r = r max = R s and since go ~ o(Z 2 ) is negligible, we get 3/2a = —1, or a = —2/3. 
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Then at y = 0, 1 ^ the condition \1/ = C is irrelevant as we said, since we can redefine 
the constant C. Then from the second (continuity) equation we get 

,2R 8 . 2 „ 2l 2 a 2 l 2 , A 

2 a + — + —- + ... =4 



' max 1 ^" 



9 r 2 18 r z 



and in the treatment of the previous subsection we had thus neglected the a 2 term, the 
equation needed to be modified, but the sign of the correction is the same. 

We can now also correct the calculation at nonzero b, by just putting the familiar cos# 
term 

o 2 b 2 I 2 a 2 

R 2 { 2p 2>{ + p 2 c 

from which we get 



% = 2* 1 + V 1 " m + m (1 + T )(1 " + -) < 5 - 59 > 



The maximum impact parameter (and thus the scattering cross section a = ^b 2 ma ^) gets also 
modified 

R 2 4/ 2 a 2 

™(i + ir)] (5-60) 



2 i? 2 

The perturbation for / ^> r (around flat 5d) will be left for future work. 



6 Aichelburg-Sexl solution in AdS background 
and scattering analysis 

In this section we will analyze the case of an A-S wave in AdS (for future application to the 
gauge-gravity duality). First, we have to derive the solution for the A-S wave inside AdS. 

6.1 Aichelburg-Sexl solution in AdS bacground 

Let us notice that [27] analyzed putting A-S Shockwaves in more general backgrounds, of 
the type 

ds 2 = 2A(u, v)dudv + g(u, v)hij(x l )dx l dx : ' (6.1) 

The calculation of the A-S solution in this background, with a source=massless photon 
at u = 0,p = was done as in flat background, just by gluing two regions at u=0 with a 
shift Av — f — f{x l ). In [27|, it was found that the Einstein equations are satisfied if 

A )V = = g >v 

A _ Qjwf = 327rp GA 2 S(p) (6.2) 
9 9 
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Indeed, in Minkowski background (A=-l/2, g=l) one finds the Aichelburg-Sexl solution, 
Af = —167rpGS^(p). Notice that if the equations are not satisfied, it just means that one 
can't find a solution for the ansatz taken. For example, spherical sourceless (p=0) waves 
of this type in flat space are excluded (A = —l/2,g = r 2 = (u — v) 2 /4 doesn't satisfy the 
conditions), but Penrose found another type of solution. 

The authors of |2Z! were able to find such Shockwaves in the Schwarzschild solution in 
Kruskal-Szekeres coordinates, 

3 

ds 2 = -32— e- r/2m dudv + r 2 (d9 2 + sin 2 9d(j> 2 ) 
r 

uv = -(r/2m-l)e r/2m (6.3) 

namely 

Notice that if one would like to put AdS in the form in (J6.1j) . one can't: For a schockwave 
moving on the brane, the AdS background would be written as 

ds 2 = — (dudv + dx\ + dz 2 ) (6.5) 
z 1 

which is not of the desired form, whereas for a wave moving in the z direction 

, o dudv + dx\ 

ds 2 = — 6.6 

(u — v) z 

which doesn't satisfy the conditions. But there could still be a solution of a different type. 
Note that neither the previous metric nor the global AdS form 

ds 2 = l 2 (-dt 2 cosh 2 p + dp 2 + sinh 2 pdttl) (6.7) 

nor the other forms 

I 2 

ds 2 = — ^{-dt 2 + d9 2 + sin 2 9dn\) (6.8) 
cos 2 9 

or (with r/l — sinhp = tan#) 

t 2 dv 2 

ds 2 = -dr 2 (l + — ) + T + r 2 dttl (6.9) 

v I 2 ' 1 + £ A v ' 

help us in putting AdS into the form desired by j2Z], so we do need something else. 

Indeed, we will see that instead we can follow closely the calculation of Emparan j3Hj, so 
we will describe it, modifying it for our purposes. 

Emparan [13] uses the metric of the one-brane RS model, perturbed with a general 
gravitational wave, in the form 

ds 2 = e~ 2lvl/l (-dudv + dx 2 + h uu (u, x\ y)du 2 ) + dy 2 (6.10) 
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But all we have to do in order to go to AdS is to replace \y\ — > y. Then, under the coordinate 
transformation 

y/l = lnz/l (6.11) 

we would get 

ds 2 = —(—dudv + dx 2 + h uu (u, x l , y)du 2 + dz 2 ) (6.12) 
z l 

which is the form that we wanted to obtain using the [27] formalism. 
But jlH] gives the Einstein tensor for the RS metric (|6.10|) as 



G - ^ d ~ l \ 
^yy ~ ;2 y w 



IX 1 

d{d-Y) 2(d- 

g^v = (— ^ 1 — S (y))9iw 

1 d„ud v u[e~ 2]yl/l (d 2 y - sgn(y)jd y ) + V 2 x }h uu (6.13) 



where we have actually corrected the jlH] result by putting the sgn(y) function. In the AdS 
case however, the sgn(y) is absent (since it came from d y \y\). 
The RS equations in the absence of h uu are 

Gab = Ag AB + X5(y)g^ B (6.14) 

(cosmological constant A in the bulk and on the brane A = brane tension) and can be seen 
to be satisfied, we could read out what A and A are. Then note that the equation for h uu is 
linear. 

In our case, adding the energy-momentum tensor of a photon of momentum p, (which 
will generate the A-S metric), travelling at fixed x % and fixed radial position in AdS, yo, 

t AB = P 6(u)S d ~ 2 (x*)5(y - y )5 A M B (6.15) 

we get an equation, with h uu = <M(-u) 

- l -[e- 2 v'\d 2 y - jd y ) + V^]$ = ^G d+lP 5 d ~ 2 {x*)5{y - y ) (6.16) 

Note that the flat space limit I —> oo gives the correct result, —l/2d 2 h = 8nGd+ipS d ~ 1 (x). 
Going to 4d Fourier space 

<S>(q,y) = J d d - 2 xe- iq - x $(x,y) (6.17) 
and similarly for t uu , one obtains 

${q, y)" - y$(g, y)' - g 2 e 2 ^$(g, y) = -lQix V G d+l 5{y - y ) (6.18) 



31 



Going back to Emparan's case [13], the previous equation would have d/l sgn(y) and 
e%l/ z . The solution to that equation in Emparan's case is 



Ae^K d/2 (e^ l lq) (6.19) 

where the Bessel function K was chosen among the 2 solutions to the Bessel equation because 
of the boundary conditions: one wanted that at y — > oo the solution dies off, not blows up 
(Id/2, the other solution, blows up exponentially at infinity). The \y\ in e d \y\l 21 was because 
of the sgn(y) in the equation, and the \y\ in the & v " 1 argument was due to the e 2 ' y '^ in the 
equation. Then both at y = oo and — oo we need the behaviour of K v (x) for x — > oo. 
Finally, the constant is fixed by normalizing the coefficient of the delta function 



A(^K d/2 (lq) + qK' d/2 (lq)) = -8nG d+lP (6.20) 



and using an identity for Bessel functions A can be put to a simpler form. Also using a more 
general energy momentum tensor for the momentum space wave, t uu (q)5(y) one has 

h U u{q,y) = MGt uu (q)e a — — - (6.21) 

qK d / 2 -x(lq) 

For the photon energy momentum tensor, going back to x space and making the angular 
integrations, using 



d-2 



Jd-4 (qr) 

dn d _ 3 e iqrcosti = tt d ^ 4 I d9sm d -*6d6e iqrcos0 = (2n)^ 2 (6.22) 

(qr) 2 



one gets 

d\y\ 



h uu (u,r,y) = - — -j^vK u )^-a dqq ? J^(qr)— — — (6.23) 

(2tt)— r~ Jo 2 K d/2 -i(lq) 

In our case, the generalization is very simple. There are no \y\ in the equation (|6.18J) . so 
none in the solution. Again the solution at y — > oo has to decay, so we choose the Bessel 
function K for y > y$. But now for yo > y — ► — oo we get the exponent of the Bessel function 
becoming K(x), x — > 0, for which K u (x) blows up as x~ v . Instead, the Bessel function 
I u (x) behaves smoothly, as x v . So the solution for y < y is with I d / 2 instead of K d / 2 . The 
normalization of the delta function is also different. 

The solution is now of the type 



$ = K d/2 (e v / l lq) y>y 

= A 2 e*I d , 2 (e yll lq) y < y (6.24) 



Continuity at y gives 



A _ hi 2 (e^ l lq) 
A 2 K d/2 (ev°/ l lq) 



(6.25) 
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and the jump in the derivative gives the delta function normalization ((A$'(?/o)) = 
— 167rGd+ipe 2ya ^). Using the Bessel function relations 

zl' v (z) + vl v (z) = zl v -x(z) 
zK' v {z) + uK v {z) = -zK v ^{z) 

h(z)K u+1 (z) + I v + x {z)K v {z) = - (6.26) 



we finally get 



h uu (q,y) = SnG d+1 t uu (q)e^e^r yo K d/2 (e^ l lq)2lI d/2 (e^ l lq) y > y 

= 87rG d+1 t uu (q)e^e^ y H d/2 (e y / l lq)2lK d/2 (e y °/ l lq) y < y (6.27) 



So now going in x space, taking the usual photon energy-momentum tensor and making the 
angular integrations we get 

h uu (u,r,y) = 8Gd V_[ p 5 ( u ) e2ie d -4 — / ^^fe(^)^d/2(e J/A /g)/d/ 2 (e ?/o// /g) y > y 



(2ir) 2 r 2 

(2< 



2 

■J U 

£rf**(«) Eg— / dqq—J^(qr)I d/2 (e^ l lq)K d/2 (e^ l lq) y < y 

2 r 2 Jo 

(6.28) 



Again, the last integration cannot be done, except on a certain hypersurface. Indeed, we 
have the relation 

/ dxx I/+1 ^(ax)/ M (te)J 1/ (cx) = — = ^V^ll (6.29) 

J M 1 M ; V 7 ^(/i 2 - l)^ 2+1 /4 v ; 

(where Q^{z) is the associated Legendre function of the second kind), that is of the desired 
form, which is however valid only if Re(a) > \Re(b)\ + \Im(c)\, Re{v) > —1, Re(fi + u) > —1 
(all satisfied) and 2ab[i = a 2 + b 2 + c 2 , which imposes a constraint. 
Thus we obtain 

(for both y < yo and y > ?/o0 on the hypersurface 

r 2 = ^eW'^e-H 1 _ i _ e — 1— ) (6.31) 

One could presumably check this by the Aichelburg-Sexl procedure, namely of boosting 
the AdS black hole and then taking the limit where the mass of the black hole goes to zero 
as the boost goes to infinity. It is however quite difficult in practice. 



h uu (u, r, y) = C-^p5{u)e-^?-«e^ C = 2 (6.30) 

(2tt)— v^tt^-I) — 



33 



6.2 Scattering analysis 

Let us look now at the AdS scattering. Let us first obtain the limits of AdS-A-S wave. 
Defining as before h uu = 3>5(u) we get 

poo 7 / 

* = C^e*&» / dzz^Ji^(z)Kw(* /l -)Id/2{^ fl -) (6-32) 

with C = 8Gd+ilp/ (2ir)^~ . As we can see, for r ^> / the integral is dominated by the region 
of small argument of I and K and we can use 

~ *f'"i>TT) ; K ^ ~ * K " {x) ~ 5 ( " " )_ " (6 ' 33) 

But 

/•oo 

dxx 2n+l J (x) = (6.34) 





so we need to expand l2{bx)K 2 (ax) up to the first term that is not of x 2n type. We find 

b 2 1 

I 2 (bx)K 2 (ax) = —^ + ct.x 2 + ct.x 4 -a 2 b 2 xHog(x) + o(x 5 ) (6.35) 

4cr 64 

and using 

poo 

J dxx 5 log(x)J (x) = —64 (6.36) 
Jo 

we get 

Instead, when r <C Z (actually, for e y l l ljr ^> 1), we can use the large argument expansion of 
I and K, 

6 X TV 

I v (x) 

\ 

and obtain (for d=4) 



I u ( x )~—=; K v {x)~J— e - x (6.38) 
\Jtv2x V 2a; 



Ce 2i l 



$ ~ ; (6.39) 

21 v / r 2 + /2 ( e y/i _ el/0 /i)2 

and so if we also have y/l,yo/l <C 1 we obtain as expected the 5d result 

C 

$ ~ — (6.40) 

2/ v / r 2 + ( y _ ?/0 )2 

Note that the result in ()6.39|) can be obtained also if r/l ~ l,e y / 1 ^> 1, which means y/Z ~ 
a few (not too large). 
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Another particular case of interest is y = y$. Then we can do the integral at all values 
of r and obtain (C = 2RJ 2 , R s is 4d the Schwarzschild radius) 



r 2 / 4Z 2 I 2 

* = 21M-1 + W2 ^""(-l + sjl + ^5^) + , ] (6.41) 



1 r 2 e -2y /l 

and we can check that for r ^> Z (and y^/l ~ 1 or <C 1) we get 

2_R S / 6 

^ — r 6 e~ 6j/() /' (6.42) 

same as the result that we obtain in this limit from the above answer for all y (7^ 7/0 )• 
We can also check that at e yo ^ l l/r > 1 we have 

e vo/i 

$ ~ IR 8 (6.43) 

r 

as we obtained from the formula at arbitrary y. 

Finally, let us now look at 't Hooft scattering in AdS$ in the two limits. For r / or 
l/r ~ 1, > 1 (so that Iq > 1 or Iq ~ 1, e y/z > 1) 

G e 2( 

$ ~ (6.44) 

21 ^/r 2 + P( e v/ 1 - e»/ ( ) 2 

and hence (since 5 = p_ and going to z = qb = qr variables and using p_p® = s/4) 

3»-»C) 

<?(6, a) = . ^ 2i 9 (6.45) 
V^ 2 + / 2 g 2 (e^ - ew/')2 

and thus if 5 is small the amplitude is 

A ~ ^^ 5Se /" — Jq(z) 



g 7 ^2 + / 2 g 2 (e^ - e»/ i ) ; 
= 2l _L e ^ fl exp[-(v / t/(e y/i - e w/ ))] (6.46) 

27T 

where the exponent is therefore large. 

However, 5 1 means either y ^ yo and G ! 4se 3 ^~ 2/0 ' ) ^ 2/ ' ) <C 1 or y ~ y and GiS-e Vo ' 1 <C 1, 
so the only possibility is y ^ yo < I, r <^ I, G^s <C 1, but we still want G^s ~ 1, but < 1 for 
't Hooft scattering, so it's not clear that there is a good regime in between. 

For r ^> I (or rather Iq <C 1), we obtain in D=4 

5(6, s) ~ 2 5 6 (6.47) 
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and therefore now 5 is always small, so 

A~—J dzz D l 2 - l J D/2 ^{z)5{z) (6.48) 
Q Jo 

Unfortunately, the result in D=4 is infinite, and to obtain the finite t-dependent piece we 
would need to get 5 at general D, which seems to be quite difficult to do, but we can say 
that the result in D = 4 + 2e will change (q/z) 6 =>- (q/z) 6+me and so 

A m = 4 ) = ^l^ e | (2y + y )^2L ln t K 0^2^ t e U2y+y ) ,g ^ 

27r2 7 3 — m 



7 Conclusions 

In this paper we have reanalyzed the question of black hole formation in the high energy 
collision of two particles via the classical scattering of two Shockwaves. 

We have found that string corrections increase the horizon area. For the effective shock- 
wave metric in [2Sj, we have found that if we scatter head-on (at b=0) two such waves, 
each characterized by an impact parameter b > R s , we obtain trapped surfaces which are 
deformed disks of area higher than the area obtained from A-S wave scattering. For the 
effective shockwave metric in j21], in the case of of R%jY 3> 1 (Y — a'log(a' s)), we get 
an increase of the area of the black hole formed, as well as of the classical scattering cross 
section, a = ^b^ iax , while in the R 2 S /Y 1 we get that the area of the formed black hole is 
of the order of Y (modified string scale), not R 2 , so much larger. 

For higher dimensions, we have found a conservative approximation scheme for the area of 
the horizon formed which gives us a maximum impact parameter (indicative of the scattering 
cross-section, as we expect that a = 7r6^ ax ). We have thus obtained that in D=4, b max = 
R s /\/2, and in D=5 for instance b max ~ 0.9523i? s , which is again a more conservative 
estimate as the one in 

What was more surprising was the fact that although graviton-graviton scattering should 
be described by the collision of two ideal sourceless waves, given in the Khan-Penrose solution, 
there doesn't seem to be a horizon forming even at zero impact parameter. There is a theorem 
that a singularity will form in the future of any sourceless wave collision, yet we can't find a 
trapped surface, namely the usual trapped surface calculation doesn't have a solution. We 
have speculated that maybe the gravitons cannot be described by sourceless waves at all, or 
maybe trapped surfaces are inherently different from the case, namely that the surfaces 
form only in the interacting region u > 0, v > 0, not at the border (u=0, v=0) as in the 
photon scattering case. 

We have extended the formalism to curved backgrounds. For more realistic scenarios, 
involving possible creation of black hole at accelerators for low fundamental scale, we have 
chosen the one brane Randall-Sundrum case. In the case that the 5th direction is highly 
curved, we have obtained just corrections to the flat 4d case, whereas for a weakly curved 
5th direction, we have corrections about the 5d flat space black hole creation. 

Finally, we have found a solution for an Aichelburg-Sexl wave inside an AdS background, 
and we have calculated the scattering amplitude for 't Hooft scattering in such a wave, at 
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small and large distances r. This was done for later use f7j for analysis of the gravity dual 
of QCD high energy scattering. 
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